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Q\ ■ Casimir interactions ( due to the massless scalar field fluctuations ) of two 

surfaces which are close to each other are studied. After a brief general 
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j> , presentation of the technique, explicit calculations are performed for specific 
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I. Introduction 

Experiments to observe and measure Casimir forces have so far been per- 
formed with the geometrical setups involving two ( actually disconnected ) 



surfaces PQ. The original parallel plate Casimir interaction is exact for in- 
finite plane surfaces j2], which in practice means that valid for planes very 
^ ■ close to each other. Effect for the parallel plane geometry were first verified in 

1958 [3]. Recently the experiment for this geometry was improved to a high 
precision |3] . The Casimir experiments other than the above mentioned ones 
have been performed for a sphere close to a plane configuration [5J: which 
do not give rise the precise alignment problem of the parallel planes. Note 
that the calculation for the sphere-plane geometry gets closed to be exact if 
the radius of the sphere is small compared to the distance to the plane [H]. 
Sphere-sphere geometry has also been studied subject to the similar approx- 
imation as the sphere-plane problem [7j. The interaction of two co-centric 
spheres has recently been addressed jS]. 

Single cavity experiments so far have not been realized [H] , which we think 
would be very interesting: For example inserting the data from quantum 
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dots ( i.e. radius ~ cm ) into the theoretical expression for the vacuum 
energy of a spherical cavity capable of confining electromagnetic field pQ, 
one gets ( in h — c — 1 units ) 0.5 ■ 10 6 cm _1 = lOeV for the Casimir energy 
which is of appreciable magnitude This is comparable to the total energy 
between the parallel plates of the latest experiment fl], i.e. E = -^-p (Area 

of Plates)~ 720(5*10-3)3 • (2 • 2)cm- 1 ~ 10 6 cm _1 . 

The purpose of the present work is to study new two surfaces geome- 
tries. We calculate the Casimir forces resulting from the vacuum fluctu- 
ations of massless scalar fields between surfaces close to each other. For 
massive fields, for any realistic experimental setups the Casimir energies are 
extremely small. The expressions always involve a factor e _AtA , where fi is the 
mass and A is the separation; which for electron and for nanometer distances 
is e ~ 2 ' 510l ° 10 7 ~ e -2500. it is practically zero. 

After a brief outline of our approach to the surface surface interactions 
in Section 2, we proceed with specific examples, that is co-axial cylinders, 
co-centric tori, co-centric spheres and co-axial conical surfaces. The case of 
co-axial cylinders my offer an experimental test in the light of the recent 
advances in stable metal nanotubes [TU] . 



II. Casimir energy for the region between two bound- 
aries which are close to each other 

We first choose the suitable spatial curvilinear coordinates rj^ , j = 1, 2, 3 for 
the geometry we deal with. The corresponding Minkowski metric and the 
Klein-Gordon operator are then 

ds 2 = dt 2 — gijdrfdrf (1) 

and ( g = det(#y) ) 

d 2 1 d f . d 



The Green function is 



Ju(\)(t-t') 



G = E 2u(x) W^wfr), ( 3 ) 



where $ w (a)(^) and w 2 (X) are the eigenfunctions and eigenvalues of the equa- 
tion for the massless scalar field 

- ^^V^J^wfo) = <" 2 (*)*«(A)fa)- (4) 



( For massive scalar field one only changes u 2 by uj 2 + /i 2 ; with being 
the mass ). We assume that the above equation is separable in the spatial 
coordinates . Here rj and A stand for the collection of the coordinates 
?7 J and the corresponding quantum numbers Xj ( which are specified by the 
boundary conditions ) respectively. The functions $oj(a)(^) are normalized 
with respect to the norm 

\\n 2 = I d 3 Vv fgm V )\ 2 , (5) 

J A 

where A is the domain of the coordinates r/- 7 . The vacuum energy density 
can be then obtained by calculating the coincidence limit derivatives as: 

T= Reg [ lim (-?- + g v -^—)G{ V , rf)). (6) 

"Reg" stands for regularization. In the specific examples it means that we 
have to subtract the terms ( in the Plana sum formulas to be employed 
over the modes ) corresponding to the vacuum energy of the free space, 
the boundary energy etc. To calculate the Casimir energy one needs the 
eigenvalues of the problem. The eigenvalues w 2 (X) depend on three quantum 
numbers Xj corresponding to the degrees of freedoms in directions in which 
we assume that the equation (j3J) can be separated. We further assume that 
after the separation of variables the eigenvalue equations in coordinates rj 1 , rj 2 
can be trivially solved, and the corresponding quantum numbers Ax,A2 are 
easily obtained. This assumption does not introduce a strong restriction. In 
fact many problems in the literature are of that type. For example when one 
studies the Casimir energy inside a spherical cavity, only nontrivial problem 
is the radial equation in which one has to deal with the roots of the Bessel 
functions to impose the boundary condition £Q. 

In this work we employ an approximation method to calculate the nontriv- 
ial spectral parameter A3, which is valid if the problem involves two bound- 
aries in direction rj 3 , which are close to each other. 
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After the separation, the problem in hand in r] 3 can be converted into the 
Schrodinger form 

<]1 +W XlX2 (r ]3 )}^ X3 (r ] 3 )=E(X)^ x . i (r ] 3 ). (7) 



L d(7] 3 ) 2 

The form of the "potential" W XlX2 (rj 3 ) and the relation between ou 2 (X) and 
E(\) depend on the choice of coordinate systems. The explicit examples are 
given in the following sections. The boundary conditions we wish to impose 
for the type of geometries under investigations are 

$A 3 (% 3 )=0, $ A3 (^)=0, (8) 

where tjq < rjf. In practice these boundary conditions require dealing with the 
roots of special functions which are quite involved. However if the boundaries 
are close to each other, instead of (jZJ) we can employ the simpler Schrodinger 
equation 

+ K*te)]*W) = (9) 

where the constant potential in the region is given by 

OO, 7] = T]q, 7] = T]'{ 

w XlX2 (^f), v'eivlvf) (10) 

The eigenvalue equation (j§J) has the following solutions 



E°(\) = (^) 2 + W° XiX2 (11) 



and 



*W) = vf 8111 ^' (12) 



where A = r/f — t]q and A 3 = 1, 2, The system given by Q is a good 

approximation if the condition 



max \W XlX2 ( V3 )-W° XlX2 \ «min|£°(A)| (13) 



is satisfied. 
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In the following sections we apply this approximation method to the spe- 
cific geometries. 



III. Casimir energy in the region between two close co- 
axial cylinders. 

In the cylindrical coordinates, i.e. with the metric 

ds 2 = dt 2 - dz 2 - dr 2 - r 2 d(f) 2 (14) 
the eigenvalue problem we have to solve is 

1 (9 <9 1 B 2 8 2 

r dr dr r 2 dcf) 2 dz 2 

After solving for the trivial coordinates z and <f) we have 

^ipz+im<j> 

$ = — — i=-v nm {r). (16) 
l~K\Jr 

Here v^ m (r) are the normalized wavefunctions corresponding to the radial 
equation 

, d 2 m 2 - 1/4 2 

t'dr 2 r 2 m ~ ^ nmVnni} \ ' 

with 



Upnm = \JV 2 + Vim- ( 18 ) 

The quantum number n should be determined from the boundary conditions 
on the co- axial cylinders with the radii r < r\\ 

v n m(r ) = 0, v nm {r x ) = 0. (19) 

The solution of (fT7j) satisfying the boundary condition at r is given in terms 
of the Bessel functions as 

Vnm{r) = ^J[i nm r , (20) 

where f2 nm to be obtained from the normalization 

n drr\v nm (r)\ 2 = 1. (21) 

ro 



In practice however the above integral is very difficult to calculate for arbi- 
trary values of r and r\. The spectrum \x nm should be determined from the 
boundary condition at n which is quite involved equation. However if the 
cylindrical surfaces are close to each other we can rely on the approxima- 
tion method summarized in the previous section. Instead of the eigenvalue 
problem (|17|) we consider the following one 

[-^ + V(r)}v° nm = (n° nm ) 2 v° nm (22) 
with the constant potential 

!oo, r = r , r = n 
J-\ °, ! (23) 

The above equation is then trivially solved as 

f2~ 

v °nm = J^sin(/i° m (r -r„)); A = r 1 - r (24) 



with the spectrum 



u° 

r^nm 



TT 2 n 2 m 2 — - 



r ri 



±; n = 1,2,3,... (25) 



For the present specific case the condition (|13|) is valid for A < ro. The 
Green function of the system is then easy to deal with: 



x x „iw° nm (t—t')+ip(z—z')+im(<l)—(l>') 



G =E E / JP - 8n2ujpnm <m{r)< m {r') (26) 



n=l m=— oo 

where 



<4nm = \P 2 + i^nJ 2 - (27) 

We insert the above Green function into the coincidence limit formula 

T = Reg [J lim {d t d t , + d r d r , + d z d z , + \d^)G], (28) 

Z t,r,z,<f>—*t',r',z',<t>' r 

where " Reg" stands for the usual regularization which will be defined explic- 
itly. The total vacuum energy per unit height is 

E= d<P rdrT = - ij; Reg^lnJ- (29) 



To perform summations we use the Plana formula pQ: 



~ . . F(0) f°° > \ f°° , F(it) - F(-it) , , 
£F(n) = -^ + / dnF(n) +i j o dt U_i ( 30 ) 



n=0 



In the application of the above formula to the summation over the radial 
quantum number n, the n = term and the second term, corresponding 
to the surface singularity and free space divergence respectively should be 
subtracted. Thus the regularization stands for (in n summation ) 



Reg E F (n)] = E F W + ^r~ dnF ( n )- ( 31 ) 

n=l n=l Z J0 



Going back to (|2^j) we first perform the sum over m. Note that since the 
argument of the square root in ()29|) is always positive we can replace it by 
the absolute value 



1 f°° dp f°° , ~ / ,7m. „ m 2 -l/4, 

E = 2 L s L rfm &9 £ v |p2 + "a + - V-i + 



00 00 



4i? 2 



2ttRE(A,^) + E 1 , (32) 



where i? 2 = rori. Here 



B ( A ^>4f Et/ll^l 2 + (^) 2 -7^l 1 133) 



''2/2' 2 7-oo (2tt) 2 L ^ V " '" ' ' v_ A y 4i? 2 
with k being 2-dim. vector k — (p, jf); and 



* - £ I E E ^ + p* + Cir- - ±\ ]]■ (34) 

It is obvious that employment of Reg J2m=i does n °t mean that there is an 
actual regularization in m summation. It simply imply the usage of (J31|) : for 
as it will be seen below that the first summation on the right hand side of 
that formula is exactly calculable. 
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To evaluate E\, we first employ (|31|) in n summation. The last term, i. 
e., the J °° dn integral term, becomes formally the same as E of (|33j) . Thus 
we can write 

E x = -2AE(nR, ^) + E 2 + E 3 . (35) 
The terms E 2 and E 3 which come from the first and second terms of (jSTjl are 

e, = i r - **9 [± \l\- 2 +v>-— 2 \ ] (36) 



and 



To evaluate £"2, we write ((Hlj) as the last term of the right hand side of (J30|) : 
then after suitable change of variables we arrive at 

1 r°° x 3 dx r°° ^ ^Jy 2 — 1 z" 1 x 3 <ir z" 00 V?/ 2 + 1 

2 ~ ~8n 2 R 2 ^Jo VTTx^Ji y e™y - 1 + 7 VT^ 2 Ji y e™y - r 

(38) 

We can easily estimate the upper limit of the above integrals. The first one 
is smaller than while the second is smaller than |. Thus 

^ < ISrW (39) 
To evaluate £3, since A < i?, we can first approximate it as: 



00 7T ^ 7i? / p 2+4# m e 27rm - 1 



2 



dp ™ dK (27rRs^p 2 + , _ 



' 7rn ^ 2 1 



2tt 2 7-oc 2tt ^ 



We use the formula [TT] 

/•OO -TT" 

/ dAAT (a;VA 2 + 6 2 ) = — e~ xb (41) 
Jo 2x 
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for the integration over p, then perform the summation over n. Finally we 
have 

which is negligible small. 

To calculate the main contribution (fHHJ) and the first term of of (135)1 we 
apply (jUI) and (jSU): 



. d 2 k r°° dn n 2 n 2 



dkk dnJ^-k 2 + fx 



»dhk r dnJ^f-k 2 + fi 2 



o 2tt Ji^/JJ^ e 2 ™ - 1 
7r 2 1 r 2 ^ y 3 rfy r 00 



1440a 3 32vr 2 a 3 



For (jHrfy i. e., with a = A and /J = since A i? we have 

1 7T 2 1 

#(A — ) (44) 

1 '2iT 1440A 3 192i? 2 A 1 ; 

Inspecting flUIJ), for a = 7ri?, /x = ^ we see that its contribution is lO™ 2 -^ 
times the second term in the above expression: thus it is also negligible. The 
final result for the Casimir energy between the close cylinders is then: 

ir 3 R M 15 A 2 
E = -T2^ (1+ 2^K ] - (45) 

Note that the inclusion of the second term in the above expression does not 
contradict our approximation of (|23jh for the contribution of the first term 
after this approximation in the potential would be of the order ^ . It is easy 
to check that in ^ — > oo limit the above result becomes same as parallel 
plate energy. 

Finally we like to remark that, for one-boundary geometries, for example 
for .D-dimensional ball there are satisfactory techniques to deal with the 
problem involving the roots of Bessel functions ^2]- We may hope that 
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these techniques may also be adopted for geometries with two-boundaries. 
For boundaries close to each other however, we can rely on the result of (|45|) . 
for it gives the correct limit of parallel plates in 4r — ► oo limit. 



IV. Casimir energy in the region between two tori 

Problem differs from the previous one by the boundary condition. Instead of 
()19|). the solution of the e- value equation (|T5)l should satisfy 

$|r=r = $|r=n = °> $ U=0 = ®\z=L (46) 

where L is the circumference of the tori. For A r we have 



$ = ^V^VA Sm( A (r - ro)) (47) 



and 



,2nk. n ,m. n ,jra,„ 1 

= V ( — > + ( ¥> + "a - «? i < 48) 



where /c, m G Z and n = 1, 2, 3, . . .; and i? 2 = r^ri as in the previous section. 
The total energy between the close tori is then, after employment of Plana 
formulas 

OO OO OO OO /*00 00 

E = - Y Re 9[J2 Y Uknm ^ = Y / dk Reg[J2 Uknm] 

m=— oo fc=— oo n=l m=— oo ^ n=l 

oo oo oo 

+ E E Reg[J2^knm) (49) 

m=— oo n=l k=l 

Note that unlike the previous case, since the degree of freedom along the tori ( 
i. e. along z-coordinate ) is also restricted, we have to perform regularizations 
for the fc-summation too. The first term on the right hand side of the above 
equation is exactly the Casimir energy for the co-axial cylinders considered 
in the previous section. Thus, we rewrite (}4~9*j) as 

E = LE C + E' (50) 
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In a fashion parallel to the evaluation of (|4U|) . taking the advantage of R 3> A 
we can evaluate RegJ^kLi i n E' : 



Id™™ , m 2 n 2 n\ , . 



From 



^ . Im 2 n 2 n 2 , r 00 , , , /m 2 7r 2 n 2 , 



A 2 



oo dmJ (2LJ%r + 



m? i 7r 2 n 2 ' 



+ 2vr / "\ v fl - — , (52) 



using the approximation e27r L_i — e 27rm for m > we arrive at [TT] 



\m A ir z n z , „.e A e 



*.(2L^ + -^J = „*(_ + ^_ ). (53) 
Thus we have 

R 9 i 1 1 1 1 
-E ~ ( 1 — ) (54) 

Since L > R ^> A this contribution is negligible small in compared to LE C . 



V. Co-axial cylindrical boxes of finite height. 

Instead of (flfi)). the solution of the e- value equation (|T5j) should satisfy 

$|r=r = $|r=n = 0, $| z=0 = $| 2 =L = (55) 

with L being the height of the cylinders. For A r we have 

sin(f z)e im * [Y~ . .nn. 
$ = A \7T sm T r - r '> 56 ) 

~K\ T'K V LA A 
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and 



^J(?£) 2 + ( ^ +( ™ )2 _J_ ; (57) 



where m G Z and n, /c = 1, 2, 3, . . .; and R 2 = r$ri as in the previous section. 
The total energy between the close cylinders is then 



OO OO OO — /*00 ^~ 

E = - Re 9[J2 X UJ knm} = ^ X / dk Reg[J2^k 

: r t i — — I , — 1 in — — ^-i m — — nv "0 , rl — 1 



OO 

/ Ah Pdl"V^ 

nml 



771— — oo fc=l m=— oo m=— oo u n=l 

1 OO OO 1 oo oo oo 

~ 7 X ^E^Onml + n X X ^E^fcfim] ( 58 ) 
m=— oo n=l m=— oo ra=l fe=l 

The first term on the right hand side of the above equality is equal to LE C , 
where E c is the Casimir energy for the co-axial cylinders geometry given 
by The third term for 4 > 1 can be explicitly calculated using its 

similarity with (f5Tj) of the previous section. Namely we have to multiply 
(JoTj) by | and make a change L — > 2L. Using (J5%j) we arrive at 

R d , 1 1 1 1 , N 

+ / . : 1755575 ) ( 59 ) 



32 dL AL _ x V?r 2 i? 2 + 4L 2 gftr^— _ , 

Since > A we can neglect the second term of the above expression. For 
L > A we have 



R ,n 1 



which is small due to the exponential factor. 

Let us consider the second term of (|58|). For the sake of simplicity we 
omit the factor — in the spectrum. Applying the Plana formula to the 
summation over m we get 

W = ^ + W, (61) 

in which 

(62) 

48R SnR 2 1 ' 

term is very small. The main contributions to the total energy than come 
from the first terms of (|58|1 and (J61|) : 

tt^L i2C(3) 
E = "72W + W (63) 
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Inspecting the above result we observe that the energy is positive around 
L < | A ( within our approximation ). Around this value of the height, the 
radial force F ra d = — |^ is repulsive. The force on the axial direction F ax i a i = 
— |f however, is repulsive for all values of L, which forces the cylinders to 
become of infinite length. When L becomes longer than |A, the radial force 
too becomes attractive. 



VI. Casimir energy between two close co-centric 
spheres. 

We employ the spherical coordinates 

ds 2 = dt 2 - dr 2 - r 2 [dQ 2 + sin 2 6d<f) 2 ) (64) 
and insert the solution in terms of the spherical harmonics 

$ = ll(M)^; Z = 0,l,2,..., -l<m<l (65) 

into the Klein-Gordon equation (0J). The resulting radial eigenvalue problem 
we have to deal is 

+ = M\n(r) (66) 

dr z r z 

subject to the boundary conditions 

vi n (r ) = 0, vi n (rx) = 0. (67) 

Here ro < r\ are the radii of the spheres and n is the radial quantum num- 
ber to be determined by the boundary conditions. To satisfy the boundary 
conditions one has to deal with the roots of the radial wave function v n i(r) 
which as in the previous section are the Bessel functions ( with m replaced 
by I + 1/2). However since we are interested in A e r] - r < r limit, we 
can proceed as we have done in the previous section. For the radial wave 
functions and the eigenvalues we obtain 



vt(r) = J — sm(ul(r-r )), (68) 
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-2 2 (1 + 1)2 

KJ 2 = ^ + li ^; n=l,2,.... (69) 

With the above approximated radial eigenfunctions and eigenvalues we can 
write the Green function as 

co oo I iu) ln (t—t') 

c = EEE -5-5— 0) 00 (70) 

n=l Z=0 m=-J 

Integrating the vacuum energy density 

T = Reg [\ lim , , fi t d t , + «9 r «9 r , + ^d e d e > + —l^-d^jG] (71) 
over the volume between two co-centric spheres we get the total energy 

CO 1 CO 

E = Y.V+~) ifeffE^- (72) 

Z=0 Z n=l 

Applying the Plana formula to the n summation and dropping the n = 
term and the integration over n we get 

E =-^mL dnF{n) (73) 

where 

00 3 

To use the Plana formula pQ 



k=0 



E/(* + 5)=y (75) 



we have to get rid off the poles of the function F(n) at the imaginary axis 
2^ns = 2mm. Thus we work with the function 

CO 3 

s ~ 2 

with /5 > 0. Then (J75J) becomes 

14 



where 

E ' = -L lim I" f SS3 I dx\lx 2 -(^) 2 ( — - + 



27rA /3-oy e 27rs + 1 V Vff V^ +is ) - 1 e x ^" is ) - 1 

(78) 

OA 

Using y <C 1 we get 

£' ~ (79) 

288A v ; 

Thus the total energy in the region between the spheres is 

7t 3 R 2 

In 2j — > oo it is obvious that the above energy approaches the parallel plate 
formula. 



VII. Casimir interactions of two close co-axial 
cones. 

The geometry we like to present in this section is two cones with common 
axis at positive ^-direction and appeces at the origin. By close cones we 
mean the appex angles 9\ and 9\ are close to each other, that is 

A = #i - 9 < ^sin #o sin 9 1 = 0. (81) 

In the above approximation the solutions we employe ( in spherical coordi- 
nates ) which vanishes at the surfaces 9 = 9q and 9 = 9\ are 

Kn = ^-J, nm ^r)^=^{ — {9 - 9 )), (82) 

where 

AW = ^) 2 + (|) 2 (83) 

The energy u in (|H2*|) is continuous. The Green function is ( with the cut off 
factor [3 ) 

oo oo I —f3u)+iu)(t—t') 

G=EE E o ®n m (r, 9, 0)$- m (r', 9', </>')■ (84) 

n=lm=0m=-l AUJ 
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Note that in this section we employ different regularization method than 
the previous ones. The cut off method is more suitable for the continuous 
energy spectra. Inserting the Green function of (jHI| into the coincidence 
limit formula and then integrating over 9 and </>, we arrive at the vacuum 
energy density at r: 

i f)2 f)2 oo oo f) 1 1 _j_ /3 2 \ 

E= -L Reg [(-^ + 2-£=) V V ^zVg±^ii±^ ] (g5) 

47rr ^L^ r2 d(32 )^ m ^ r 

-Reg^ stands for the cut off regularization, that is we pick the finite part of 
the expression in f3 — > limit. In deriving we used the formula [TT] 

/ ^e-^(J.K)) 2 = — Q,-i/ 2 (l + ^), (86) 

where Q u (x) is Legendre function of the second kind. We rewrite the expres- 
sion (|83|) as 

1 oo oo 

E= ^%[OE E Q-i/ 2+ ^ m (l + 2/ 2 )] (87) 

n=l m=-oo 

with 

~ (9 y 2 + 1 <9 2 

Applying the Plana formula to the summation over m we arrive at 

E = E + E U (89) 

where 

I ^ °° ^00 

£ = — ite(/ tf [0 2 / rfmg_ 1/2+/inm (l + y 2 )] (90) 



and 



1 ^ 00 oo rfmtanh J(f ) 2 - (^ 2 ) 

- Jfe^O E y. en ^b^l P ^^v^f)MW (1 + * )] 

' n=l A o x v 

(91) 



Making use of ZL ^ 3: ^> 1, tanhx < 1 and 



Reg y [dP^ /2+ls (l + y 2 )] = l-j (92) 
16 



we get 

1 00 roo 7 ( — ) 2 - ( — ) 2 

W * $ E /.e. - — 2 — )l- ( 93 > 

That is 

l^il < -^e- 2 -^. (94) 

' ' ~ 47rAr 4 v ; 

Thus Ei is negligible small. To evaluate Eq, we apply the Plana formula to 
the summation over n in (|90|). The formula we obtain is 

B ° = £ + ^-^ (95) 

where 

9A p /-°° tanhyy-z* 7 2 2 

A = Reg y [0 / dssQ_i +s (l + y 2 )] (97) 

^ POO 

B = Reg y [0 dsQi +s (l + y 2 )] (98) 
J 

Changing the variables y = t, x 2 = t 2 — k 2 ()96jl can be rewritten as 

0A r 00 dtt r 1 dkk , w 7 , 9 9n . . 

£ = -/ —— / tanh fct --fc 2 t 2 99 

Inspecting the integrals over k, that is the terms 

fi(t)= r i^ tanh(fet) (100) 
jo VI - ft 

and 

M*) = f 1 -}St=ta,nh(kt) (101) 
Jo vl — ft 

we see that both approach very fast from the value /i(0) = /2(0) = to their 
respective asymptotic values fi(t — > 00) = 1 and fzit — > 00) = 0.66 = |. 
Let us treat the second term in in detail. We approximate / 2 (i) as 

/.(«) = { f t ^ (102) 
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where a and b are both of order 1. The second term in (|99|) then becomes 



6A . r b dtt 4 f°° dtt 3 , 
^ = -7T=3(a /_ ^ + / ^) ( 103 ) 



2vrr 4 v Jo e 2At - 1 A e 2A * - 1 

Since A < 1, we can approximate the denominator of the first integrand 
as e 2A * — 1 ~ 2Ai. In the second integral making the change of variables 
2At = s , we can replace the lower boundary as 26A ~ 0. Thus ()103|) 
becomes 

_ QA ab 4 1 y°° dss 3 Qab 4 Btt 3 

2 ~ ~27rr 4 te + 16 A 4 Jo e s - V ~ 167rr 4 720r 4 A 3 ' ^ ^ 

It is obvious that the first term is negligible compared to the second. Similar 
treatment shows that that the first term in (}9"9"j) gives contributions of orders 
0(A) and O(^) both are small. Inspecting (JDTij) we see that the second and 
third terms in E are also negligible. Thus the final result for our Casimir 
energy is : 

B ^-725^ (105) 

Note that the above "density" is an expression obtained after integrating 
over G and 0. If we divide (|1(J5|) to the angular integral 



/■27T 

sin# / d(j) ~ 27t6A (106) 

1 6 JO 

we obtain the energy density averaged over the angular variables: 

7T 2 

E~£ 2 = — - + 0(A- 3 ). (107) 

1440r 4 A 4 v ; K ' 

In small A limit the above result is in perfect agreement with the energy 
density in the region between two infinite planes with angle A between them 
( i. e. the wedge problem ) jTSj 

^ = "il4o7w ( ^~^ ) - (108) 
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